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High Mass Star Formation by Gravitational 
Collapse of Massive Cores 

By MARK R. KRUMHOLZH 

^ Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544-1001 

The current generation of millimeter interferometers have revealed a population of compact 
(r < 0.1 pc), massive (M ~ 100 Mq) gas cores that are the likely progenitors of massive 
' stars. I review models for the evolution of these objects from the observed massive core phase 

I through collapse and into massive star formation, with particular attention to the least well- 

understood aspects of the problem: fragmentation during collapse, interactions of newborn stars 
with the gas outside their parent core, and the effects of radiation pressure feedback. Through 
a combination of observation, analytic argument, and numerical simulation, I develop a model 
QQ ■ for massive star formation by gravitational collapse in which massive cores collapse to produce 

single stars or (more commonly) small-multiple systems, and these stars do not gain significant 
mass from outside their parent core by accretion of either gas or other stars. Collapse is only 
very slightly inhibited by feedback from the massive star, thanks to beaming of the radiation by 
^ ' a combination of protostellar outflows and radiation-hydrodynamic instabilities. Based on these 

findings, I argue that many of the observed properties of young star clusters can be understood 

(N . 

as direct translations of the properties of their gas phase progenitors. Finally, I discuss unsolved 
I problems in the theory of massive star formation, and directions for future work on them. 

o 

. 1. Introduction 

^ ■ 

Q^. Massive star formation occurs in the densest, darkest parts of molecular clouds. These 

clumps of gas have masses of thousands of Mq, radii <J 1 pc, volume densities of ^ 10^ H 
atoms cm"'^, column densities of ~ 1 g cm^^, visual extinctions of hundreds, and velocity 
^ I dispersions of several km s~^. Observations often reveal indicators of massive star for- 

■ mation such as maser emission and infrared point sources within them, but the majority 

of their mass appears to be dark and cold. Due to their high extinctions and low tem- 
peratures, these regions are only accessible to observation through millimeter emission , 
Itn] I either in molecular lines (e.g. Plum e et '^119 97': 'Shirlev et al'2003; 'Yonekura et a;."2005t 

IPillai et al. 2006 ftl or dust conti nuum (e.g. .Carey et al, 2000: Mueller et al. 2002.1. or 
through infrared absorp t ion (e.g. lEgan a/.l]l998l: Menten et an 2005HR,athborne et a,l\ 
1200,'THSimon et oiJEoOft lR,athborne et ai!.ll200fil) ^hev are hkely the progenitors of the 
rich clusters in which massive stars form. 

In the last few years observations using millimeter interferometers to obtain still higher 
resolution have identified "cores" within these dense clumps, objects small enough that 
they approach the stellar mass scale. Cores are distinguished by even higher volume 
densities than the massive clumps clumps around them, 10^ H cm^'^ or more, and smaller 
radii, r ^0.1 pc. Some show mid-IR point sources in their centers (e.g. IPillai et a,l\ 
12006 al) . while other show no MIR emission or even MIR absorption, indicating that they 
are starless or contain only very low mass stars (e.g. ISridharan et aZ.1 120051) . In some 
cases they show signs of molecular outflows but not MIR emission, indicating that the 
extremely mass ive core contains a ve ry low mass protostar, and thus is near the onset of 
star formation l|Beuther et a^JlioOSfcll . 

The characteristic mass, size, and density of massive cores make them appealing can- 
didates to be the progenitors of massive stars (e.g. .Garavii2005j) . Moreover, as I discuss 
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in more detail in § [3 the observed mass and spatial distribution of protostellar cores is 
quite similar to that of stars in young clusters. If massive cores are the direct precursors 
of massive stars, then one can explain many of the properties of newborn clusters directly 
from the observed properties of their gas phase progenitors. The goal of this review is 
to construct a rough scenario based on this idea by following observed cores through 
collapse, fragmentation, accretion, and feedback, to the final formation of massive stars. 
In § 121 1 briefly review observations of the properties of massive cores to provide initial 
conditions for this scenario. In § |21 §0] and §1211 discuss three major questions about 
how these cores turn into stars: do they fragment into many stars or only a handful? Do 
the stars they form accrete significant mass from outside the parent core? Does feedback 
significantly inhibit accretion? Finally, in § (HI I discuss some outstanding problems in the 
modeling of massive core evolution, and suggest directions for future work. 



2. Massive Cores: Initial Conditions for Massive Star Formation 

We know disappointingly little about massive cores, despite great observational effort. 
Due to their small sizes and large distances, massive cores are only marginally resolved 
even in observations with the highest resolution telescopes available. Nonetheless, obser- 
vations do allow us to determine some gross properties of individual massive cores, and 
of the massive core population as a whole. Observations show that cores are centrally 
concentrated, although the exact density profile is difficult to determine with interferom- 
eter measurements, and fairly round, with aspect ratios of roughly 2:1 or less. They are 
cold, T w 10 — 40 K, except near stellar sources, so their observed velocity dispersions of 
^ 1 km imply th e presence of highly supersonic motions (e.g. iReid fc Wilsonll2005t 
iReiither et >^l2005ftt [2'o 06). At the characteristic density of ~ 10^ H cm ^ found in these 
cores, the free-fall time is only ~ 10^ yr, so the implied accretion rate when a massive 
core collapses is 10""* — 10"'^ Mq yr~^. 

iMcKee fc T an (,200.^ propose a simple self-similar model of massive cores in which the 
core density and velocity dispersion are power law functions of radius such that at every 
radius turbulent motions provide enough ram pressure to marginally support the core 
against collapse. The central idea is that, at the high pressures found in massive star- 
forming regions massive cores must be supported by internal turbulent motions. While 
a self-similar spherical model is obviously a significant simplification of a turbulently- 
supported gas cloud, it provides a reasonably good fit to the available observations, and 
makes it possible to calculate quantities such as the time scale for star formation and 
the relationship between core mass, column density, pressure, and velocity dispersion. It 
also provides a good starting point for simulations and more detailed analytic work. 

For the massive core population as a whole, we know somewhat more, and the obser- 
vations bolster the idea that massive cores might really be the progenitors of massive 
stars. Several authors using different techniques and observing different regions find that 
the mass distribution of massive cores matches the stellar initial mass function, shifted 
to higher mass by a factor of a few, with a Salpeter slope of F « —2.3 at high masses and 
a flattening at lower masses (Beuthcr fc Schilkc 2004: Reid fc Wilson 2005, 2006a 
This extends earlier observational work indicating that in low and interr nediate mass 
star-form ing regions the core mass function resembles the stellar IMF (e.g. iMotte et all 
ll998tlTtes ti & Sargent 1998; Johnstone et al. 2001; Onishi et al. 2002), and suggests that 
the stellar IMF may be set simply by the mass distribution of prest ellar cores, reduced 
by a factor of a few due to mass ejection by protostellar outflows ijMatzner fc McKe3 
12000.') . Simulations and analytic arguments, can, in turn, explain the core mass distribu- 
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tion as arising na turally from the supersonic turbulence present in star-forming clumps 
( Padoan fc Nordlund 2002 : Tillev fc Pudritz...2004: Li et al. 2QQ4) . 

IClark fc Bonnelll 1)200611 argue that the mass distribution of bound cores in simula- 
tions does not have a Salpeter slope and thus is unlikely to be the origin of the stellar 
IMF. However, this misses a critical point: the Salpeter slope is only observed for stars 
significantly ab ove the peak of the IMF. The f ull IMF is closer to a broken power law 
l|Kroupall200l or a lognormal l)Chabrieil 12003') . with the break or peak at ^ 0.5 Mq. 
This is roughly the Je ans mass in star- forming clumps, and indeed the simulations of 
IClark fc Bonnelll l)2006|) do show something like a lognormal distribution, with a peak at 
roughly the Jeans mass of their simulations. (The simulations are scale-free, since they 
include only hydrodynamics and gravity and use an isothermal equation of state.) 

Furthermore, recent observations fo cusing on the spat ial distribution of cores have 
shown that cores are mass segregated (|Stanke et QiJ l20o3) in much the same manner as 
stars in very young clusters: the core mass function has the same lognormal or broken 
power law form everywhere in clump, with the exception that the most massive cores, 
those larger than a few Mq in size, are found onl y in the ve ry center. The stellar popula- 
tion of the ONC follows the same pattern (Hillen brand fc Hartmann.l998.:,Huff fc Stahleil 
indicating that the observed mass segregation in stars may simply be an imprint- 
ing of the prestellar core spatial distribution. At least some of the mass segregation must 
be primordial rather than a result of dynamical evolution (i,Bonnell fc D avies 1998), al- 
thoug h recent evidence that cluster formation takes several crossing times l|Tan et all 
12006') suggests that evolution may be important too. Nonetheless, it is quite suggestive 
that both the IMF and the spatial distribution of stars in a cluster seem to be explicable 
solely from the observed distribution of gas from which star clusters form. However, the 
origin of the mass segregation of cores is at present unknown. 



3. Fragmentation of Massive Cores 

It is only possible to explain the properties of stars in terms of the properties of cores if 
there is a more or less direct mapping from core mass to star mass. Such a mapping exists 
only if cores do not fragment too strongly, i.e. if massive cores typically produce a one 
or a few massive stars, rather than many low stars. Fragmentation to a few objects does 
not present a problem, since observationally-constructed mass functions are generally 
uncorrected for multiple systems, but fragmentation to many objects does. 

One might expect massive cores to fragment because they contain many thermal Jeans 
masses of gas. At the densities of 10^ H cm"'^ and temperatures of ~ 10 K typical 
of massive cores, the Jeans mass is only Mq, so one might expect a 50 M© core 
to form tens of stars. Dobbs et al. (2005) simula te the collapse and f ragmentation of 
massive cores with initial conditions based on the iMcKee fc TanI (|2003>) model, using a 
code that includes hydrodynamics and gravity. They try both isothermal and barotropic 
equations of state. (Barotropic here means that the gas is assumed to be isothermal at 
densities below some a critical density, chosen to be 10^ g cm^'^ in the Dobbs et al. 
simulations, and adiabatic at higher densities). Dobbs et al. find that the cores fragment, 
forming anywhere from a few to several tens of objects, depending on the assumed initial 
conditions and equation of s t ate. I n no case do their simulations from a massive star. 

However, the lDobbs et all ( 2005|) calculation omits the critical effect of radiation feed- 
back from the forming star. The high densities in massive cores produce high accretion 
rates, so that the first protostar to condense within a core will immediately produce a 
large accretion luminosity as the gas that falls onto it radiates away its potential energy. 
For a typical accretion rate of lO"** Mq yr~^ at massive core densities, a 0.1 Mq, 1 Rq 
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Figure 1. Temperature {left panel) and Jeans mass {right panel) versus radius in a core with 
a mass of 50 Mp), a column densit y of 1 g cm~'^, and a density profile p oc r~^'^, taken from the 
models of lKrumholz et al\ i2006afl . The lines show the result from a radiative transfer calculation 
when the central protostar is 0.05 Mq {thin solid line) and when it is 0.8 Mq {thick solid line), 
and from using the barotropic approximation {dashed line) or an isothermal equation of state 
{dotted line). The Jeans mass is computed using the density and temperature at each radius, 
and is defined as Mj = ^■K^^^[kBT/{Gfi)]^^^ p~^^^ , where p is the density, T is the temperature, 
and fi = 2.33mH is the mass per particle for a gas of molecular hydrogen and helium in the 
standard interstellar abundance. 



star releases approximately 300 Lq of accretion power. Because the core is very optically 
thick, the radiation is trapped within it and heats the gas as it diffuses out. As a result, 
the densest, inner parts of the core where fragmentation is most likely to take place 
are subject to rapid heating, which suppresses fragmentation. Isothermal or barotropic 
a pproximations completely ignore this effect. 

iKrumholz et a,l\ 1 2006^ ) examines how feedback heating affects fragmentation in the 
context of a simple mo del of core accretion u sing a high accuracy analytic radiative 
transfer approximation llChakrabarti fc McKedEool . Figure Q s hows a sample result , 
the radial temperature profile and Jeans mass versus radius for a iMcKee fc Ta^ ll200,'^ 
core with a mass of 50 Mq and a column density of 1 g cm~^ accreting onto a protostar 
in its center. The figure compares the results using a radiative transfer approach to 
what one would find by neglecting radiative transfer and simply using a barotropic or 
isothermal equation of state. As the plot shows, both an isothermal equation of state and 
the barotropic approximation make order of magnitude errors in the temperature and 
Jeans mass. 

One might worry whether this analytic treatment done in spheri cal symmetry applie s 
to more realistic massive cores with complex density structures fe.g. lBonnell et al\\200(^ . 
The natural way to address this question is with radiation-hydrodynamic simulations of 
massive core evolution. Krumholz, McKee, & Klein (2006, in preparation) simulate the 
collapse and fragmentation of massive cores to examine this effect. The simulations use 
an adaptive mesh refinement radiation code to solve the Euler equations of gas dynam- 
ics coupled to gray radia tion transport and ra diation press ure force in the flux limite d 
diffusion approximation ^Tnielove et aliflflflsl (Kleinlll flflflt iHoweH fc Green on ghll2o"oi. 
They use the adaptive mesh capabi lity to guaran tee that the local Jeans length is always 
resolved by at least 8 cells ( Trucl ove et a/jll997l) . and that the radiation energy density 
changes by no more than 25% per cell, so radiat ion gradients are well resolved. The code 
uses Eulerian sink particles to represent stars (|Krumholz et giJl200 4h. and the sink par- 
ticles are in turn coupled to a simple protostellar evolution model ()McKee fc Tan.200,'^) 
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which computes the instantaneous stellar luminosity, including the effects of accretion, 
Kelvin-Helmholtz contraction, deuterium burning, and hydrogen burning. This luminos- 
ity becomes_a^omce_tem the radiation equation. Further details on the code are given 
in lKrumholz et all ||2005qI) . 

The simulations begin with cores following the model of iMcKee fc TanI l|2003l) . The 
initial density profile is chosen with p cx r^^^^, to a maximum density of p = 10^^* g 
cm~^, corresponding roughly to the density of the inner, thermally-supported zone of 
McKee & Tan cores. The temperature is 20 K throughout the c ore. There are initia l 
turbulent velocities chosen from a Gaussian random distribution ijPubinski et 
with a power spectrum P{k) oc k~'^(Pk over wavelengths ranging from the size of the 
core to the size of the inner thermal zone, subject to the constraint that the initial 
velocity field be divergence-free. The magnitude of the velocity fiel d is normalized to giy e 
approximate hydrostatic balance on the largest scale (equation 6 of lMcKee fc Tanll20n,'^ . 
The simulations reach a maximum resolution of 10 AU. 

Figure |2] shows the column density distribution in a simulation of a core with an initial 
mass of 100 Mq and radius of 0.1 pc, 2.0 x lO'' yr (0.37 mean-density free-fall times) 
after the start of a simulation. The core is not forming many stars, it is forming a triple 
system. Moreover, it is a highly unequal triple: the masses of the three stars are 5.33 
Mq, 0.31 Mq, and 0.16 Mq, so the vast majority of the mass has gone into the most 
massive object, the one at the center of the large disk. There are no apparent signs of 
further fragmentation, so unless feedback disrupts this system, it seems destined to form 
a massive star incorporating a significant fraction of the initial core mass rather than 
dozens of small stars. The weak fragmentation wc find from simulations provides strong 
support for the idea that the core mass function directly sets the stellar mass function. 

To understand the origin of the weak fragmentation, it is helpful to examine the tem- 
perature distribution in the core. Figure |31 shows the temperature distribution in the 
simulation dX t — 1.6 x lO"' yr, when the central star is only 3.83 Mq- At this point the 
star has not yet begun hydrogen burning, and the luminosity of a few thousand is 
entirely due to accretion. This accretion power has doubled the initial temperature of 
the gas out to more than 2000 AU from the central star, and increased the temperature 
to more than 100 K over a radius of many hundreds of AU. This heating strongly sup- 
presses fragmentation in the densest gas, where it is most likely to occur. Of the two 
stars that do form in addition to the most massive, one does so at an initial separation 
of several thousand AU, far enough that it can condense, and the other does so inside 
the protostellar disk, where the high column density provides shielding against the stellar 
radiation and produces a lower temperature. In examining a movie of the simulation, one 
clearly sees many overdense clumps that look like they might collapse, but do not do so 
because they are bathed in the radiation field of the central star. Rather than forming 
stars of their own, they fall onto the central star and accrete. 

Part of the reason suppression of fragmentation is effective is because, despite the 
complex density structure shown in Figure |21 the temperature distribution is relatively 
round and smooth. The only significant deviation from sphericity is in the protostellar 
disk. Thus, the gas around the protostar is heated quite uniformly, and outside the 
disk, where shear suppresses most fragmentation, there are no cold spots favorable to 
fragmentation. This is to be expected: the entire core is very optically thick, so radiation 
diffuses outward rather than free-streaming. As a result, there is very little shadowing, 
and clumps that are only starting to collapse are not sufficiently overdense to exclude the 
radiation field and remain cooler than their surroundings. As I discuss in §[S1 only when 
gas reaches the densities typical of accretion disks or when optically thick structures 
begin to form can there can be significant temperature anisotropics due to coUimation. 
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Figure 2. Column density of the entire core {upper panel) and zoomed i n on the forming 
stars {lower panel) for a simulation of a 100 Mq, 0.1 pc lMcKee fc TanI |2003) core at a time of 
2.0 X lO'' yr. The positions of the stars are indicated by the diamonds. Their masses are, from 
left to right, 0.31 Mq, 5.33 Mq, and 0.16 Mq. 

The weak fragmentation shown in simulations with radiation is strikingly different from 
what one obtains without it, where the number of fragments generally approaches the 
number of thermal Jeans masses in the initial cloud. Figure 21 shows why: the barotropic 
equation of state severely underestimates the temperature, making fragmentation far 
easier than it should be. The magnitude of the underestimate ranges from factors of 
a few at distances of thousands of AU to orders of magnitude in the central hundreds 
of AU. Since the Jeans mass depends on temperature to the 1.5 power, the error in 
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Figure 3. Tempera ture distribution in a 2D slice through a simulation of a 100 Mq, 0.1 pc 
[McKcc fc TanI (120031) core at a time of 1.6 x 10* yr. The image is centered on the 3.83 Mq 
star, indicated by the diamond, the most massive present in the simulation at that time. The 
outermost contour corresponds to a temperature of 40 K, and each subsequent contour represents 
a factor of 2 increase in the temperature. 



the critical mass for fragment growth is larger still. It is easy to understand intuitively 
why the barotropic approximation fails so badly: the physical assumption underlying 
the barotropic approximation is that above some critical density the gas cannot radiate 
efficiently, and all its gravitational potential energy is converted into heat. However, 3D 
simulations of star formation cannot resolve stellar surfaces, so any gas that falls into 
sink particles of radius ~ 10 AU disappears from the simulation, taking its gravitational 
potential energy with it. However, since potential energy varies as the vast majority 
of the energy is released in the final plunge from ~ 10 AU to the stellar surface. In the 
barotropic approximation one simply ignores this energy, which is the dominant source 
of heating until nuclear burning begins. 



4. Competitive Accretion 

Weak fragmentation means that a significant fraction of the mass in a massive core 
will end up in a single star or a few stars. However, for cores to be the direct progenitors 
of massive stars, it must also be the case that any additional mass a star accretes from 
outside its parent core is small compared to the stellar mass. The idea that most of a 
star's mass comes not from a parent core, but from gas in the cluster-forming clump 
that was not originally bound to that star, is called competitive accretion (jBonnell et a]] 
l200fi[ and references therein) . A number of numerical simulations appear to show exactly 
this process, and several authors have made simple theoretical models based on these. In 
these models, all stars are born from cores at roughly equal masses, with the initial mass 
ranging from brown dwarf masses ijBate fc B onncU 2005 ) to as rn uch as ^ 0.5 M©, the 
peak of the IMF, in the most recent models (fBonnell fc Batell200fil) . In these most recent 
models, most of these seeds do not accrete much mass in addition to that in their parent 
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Figure 4. Ratio of the temperature distribution sh own in Figure El to the temperature one 
would infer using the barotropic equation of state of lDobb~e t al. (200^. Note that the color 
scale starts at a ratio of 1.0, so any color other than black indicates that the true temperature 
is higher than the barotropic temperature. The outermost contour corresponds to a ratio of 0.4 
dex, increasing by 0.2 dex with each subsequent contour. 



core, but a few stars are in special locations that allow them to undergo rapid accretion, 
reaching high masses. This process determines the IMF above the peak. 

4.1. Under What Conditions Does Competitive Accretion Occur? 

Competitive accretion definitely occurs in some simulations, and there is no reason to 
doubt that those simulations produce the correct result for the physics they include. 
However, in order to determine whether the simulations accurately reflect the properties 
of real star-formng clumps, it is necessary to investigate the physics behind the com- 
petitive accretion process. Krumholz et al. (2005(i) defines the fractional mass change 
/m = A/»tdyn/A/* as the fractional mass change that a star of mass undergoes per 
dynamical (crossing) time ^dyn of its parent clump, where M» is understood to refer to the 
accretion rate after the star has consumed its initial bound core. Competitive accretion 
models require /m S> 1. Accretion of gas that is not initially bound to a star can occur 
in one of two forms: either the star may capture other gravitationally bound cores and 
then accrete them, as proposed for example by Stabler et al. (20Qfl}, or it may accrete 
gas that is not organized into bound structures. 

The former process is reasonably easy to understand, since it is simply an extension of 
standard calculations of collision rates in stellar dynamics. The only significant complica- 
tion is that collisions between stars and cores that occur at relative velocities that are too 
large do not result in capture, since the star will simply plough through the core without 
dissipating enough energy for the two to becom e bound. Even with thi s complication, 
the calculation is relatively straightforward, and lKrumholz et'al\ l(2005dl) show that the 
fractional mass change due to captures of cores with mass comparable to the stellar mass 
M* in a star-forming clump of mass M is 



/m- 



cap 



O.40CO [4 + 2u^ -{4 + 7.32u2) exp{-1.33u^)] 



(4.1) 
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where 4>co is the fraction of the clump mass that is in bound cores, u ~ 
is the ratio of the escape velocity from the surface of a core to the velocity dispersion 
in the clump, and avir is the virial parameter for the core, roughly its ratio of turbulent 
kinetic energy to its gravitational potential energy. The significant thing to notice about 
this expression is that it does not approach unity unless u is quite large, which in turn 
only happens for virial parameters Ovir <C 1, i.e. for clumps where the turbulent velocity 
dispersion is small compared to what is needed to prevent collapse. 

Accretion of unbound gas is somewhat more complex, since to determine the accretion 
rate one require s a theoretical model for Bondi-Hoyle accretion in a turbulent medium. 
iKriimholz et all ( 2005 fei 12006 have developed such a theory and shown that it repro- 



duces the results of simulations quite well. Figure shows a samp le of an adaptive mesh 
refinement simulation in which a grid of 64 Eulerian sink particles (jKrumholz et a?j|2004fl 
are placed into a turbulent medium and allowed to accrete until the mean accretion rate 
reaches equilibrium. Figure El shows a comparison of the model prediction for the prob- 
ability distribution of accretion rates to the simulation results. The model predicts, and 
the simulation confirms, that the mean accretion rate for a star of mass accreting 
from a medium of mean density p and ID velocity dispersion a is 

M,4^0BH-p ; ^ , (4.2) 

where the quantity 4>bh is a function of the Mach number and size scale of the turb ulent 
region, an approximate analytic form for which is given in lKrumholz et "dl l)20066l) . For 
the properties of observed star-forming regions, it is generally <^ 5. From this result, one 
can compute /m due to accretion of unbound gas in a star-forming clump of mass M: 

/m-BH « lO0BH<i?(A^*A^)- (4.3) 

Again, for cluster-clumps hundreds to thousands of Mq in mass, /j^j^bh can be of order 
unity only for avir *C 1- 

Combining the two potential sources of mass, one can derive an approximate criterion 
that a star-forming gas clump must satisfy in order for competitive accretion to occur 
within it. For seed stars of mass 0.5 Mq, this condition is 



,2 



Af^SOMg. (4.4) 



Straightforward application to observed star-forming clumps shows that they are nowhere 
near meeting this condition, since their typical masses are many hundreds to thou- 
sands of Mid, and their observed virial parameters are generally near unity. From this, 
^rumholz et al. ( 2QQbd^ conclude that competitive accretion does not occur in real 
clumps. It occurs in simulations only because in those simulations there either was ver y 
little turbulence present in the initial conditio ns Ce.g. Ixiessen fc Burkertll200d l200l|) . 
or the initial turbulence has decayed away fe.g. iBonnell et qiJl200,'^) . leaving the clumps 
sub- virial. 

In light of this claim, iBonnell &: Batd l)2006(l have re-analyzed the simulations of 
^onnel l et al. ( 2003) . They argue that competitive accretion occurs for a few stars that 
sit at the center of collapsing regions. These regions are stagnation points of the larger 
turbulent flow, where the velocity dispersion is much smaller than the mean and the den- 
sity much larger. Thus, although the clump as a whole is turbulent, the regions where 
massive stars form are effectively decoupled from the large-scale flow. Since these decou- 
pled regions have masses much smaller than that of the entire clump, and their virial 
parameters are smaller than unity due to their state of global collapse, they satisfy the 
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Figure 5. Slices thro ugh a simulation of Bondi-Hoyle accretion in a turbulent medium by 
iKrumholz et al.{ i2006fef) . showing density in the entire simulation box {upper panel) and in a 
small region (indicated by the black box) around one of the accreting particles (lower panel). 
The position of the particle is indicated by the small white diamond. The density and length 
are in dimensionless units where the mean density in the box is unity, and the box extends from 
— 1 to 1. The maximum resolution of the simulation is 8192^. 

competitive accretion condition H4.4(l . and stars within them gain mass via competitive 
accretion. 

I discuss whether such decoupled regions exist in real clumps below, in 14.21 However, 
one potential difficulty with the idea that competitive accretion occurs in locally detached 
parts of the flow is that it assumes that such detached regions would fragment and still 
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Figure 6. Cumulative {left panel) and differential {right panel) distribution of accretion rates 
measured for particles in the simulation shown in Figure |S] Accretion rates are normalized to 
Mo = 47rp(GMt)^/ {y/Sa)^, where p is the mean density, M, is the stellar mass, and a is the ID 
velocity dispersion in the simulation box. The histogram shows the simulation results, the solid 
line shows the Krumholz et al. ( 2006 M model, and the dashed and dotted lines show alternative 
models. 



produce many stars to compete, as it does in non-radiative simulations. However, these 
collapsing regions bear a striking resemblance to the objects we observe as massive cores: 
they are bound, high density regions with coherent velocity structures. As shown is 
jjEl such objects do n ot fragment strongly when one includes radiative transfer. Indeed, 
iBate fc BonTeiil l|2005|) estimate the typical size scale of these collapsing regions as ~ 1000 
AU, well inside the effective heating radius of a single rapidly accreting star. Instead of 
competitive accretion, the true behavior might simply be monolithic collapse of a massive 
core to a single system. In this case competitive accretion would be nothing more than 
core accretion simulated without a sufficiently accurate treatment of radiative transfer. 
Investigating that possibility will have to wait for future simulations. 

4.2. Competitive Accretion and Clobal Collapse 

ThelBqnnell fc Batel l)2nO(T|) analysis of competitive accretion concurs with lKrumholz et a,l\ 
\ 2005 J ) that competitive accretion cannot occur under typical conditions in a star form- 
ing clump. Instead, the possibility of competitive accretion depends on the existence 
of long-lasting stagnation points within which the gas and stars move together with 
gravity as the only significant force. The existence or non-existence of such stagnation 
points is something that is at least indirectly subject to observational determination. As 
iBonnell fc Bate (2006) point out, if any source adds a significant amount of turbulent 
kinetic energy to the star-forming clump, such that the rate of energy injection is com- 
parable to the rate at which energy is lost through radiative shocks, then ram pressure 
from flowing gas will push on gas but not on stars. No long-lived stagnation points will 
exist, and gas and star velocities will become decoupled. In this case, stars will ran- 
domly sample the density and velocity field of the clump, consistent with the treatment 
of Krumholz et al. (2005rf) and Krumholz et al. (2006&), and competitive accretion will 
not occur. 

The most likely candidate for significant energy injection on the size scales of cluster- 
forming molecular clumps is feedback from protostellar outflows (,Norman fc Silk.l980(l . 
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Observational efforts to estimate wlietlier tlie kinetic energy added by outflows is sig- 


nificant tentatively conclude that it is (Williams et ai. 




.Guillen et al. 


2005). While 



these results are preliminary, if they are confirmed then competitive accretion cannot 
occur. One can also look for signs of large-scale collap se onto stag nation point s. While 
there are a few examples of apparent infall signatures ijMotte et al. 200 5: P eretto et all 
on the scale of thousand Mq star-forming clumps (as opposed to onto individual 
protostars, which is expected whether competitive accretion occurs or not), the major- 
ity of searches for infall signatures have turned up negative ( G.aray 2005) . Competitive 
accretion predicts that infall should be ubiquitous. It is, however, possible that these 
non-detections are due to observational confusion or lack of resolution. 



A potentially more promising approach is to look for indirect signs of significant energy 
injection from feed back. If not, then the turbulence in protocluster gas clumps should 
decay rapidly fe.g. IStone et a/.lll998(l . leading to global collapse in roughly a crossing 
time. Over this time, a substantial fraction of the gas should be converted into stars. We 
can therefore test whether competitive accretion can occur by observationally estimating 
the time scale of star cluster formation and the star formation rate in pre-cluster gas. 

On the question of time scale, lElmegreenl l)2000l) has argued that star formation pro- 
ceeds in a crossi ng time, but a more careful examination of much of the same data by 
IXan et al\ l)2006() shows that the typical time scale is closer to 3 — 5 crossing times (6 — 10 
free-fall times). There are several lines of evidence for this, including the fairly round 
morphologies of most protocluster clumps, the lack of strong sub-clustering apparent in 
young clusters such as the Orion Nebula Cluster, the estimated momentum carried by the 
combined protostellar outflows of forming clusters, the age spreads of stars determined 
by fitting to pre-main sequence evolut ionary tracks, a,nd th e ages of dynamical ejection 
events. In contrast, in the simulation of lBonnell et o,l where feedback is neglected, 

star formation has completely ended and more than 50% of the gas has been accreted in 
1 crossing time, ^ 5 x 10^ yr. Figure shows an example of how this prediction for the 
time scale of star cluster formation compares to observations: the solid line shows the 
observed ages of stars in IC348, a cluster whose most massive member is a mid-B star, 
computed by Palla & Stable r (2000). The da shed vertical line shows the crossing time 
of « 0.6 Myr, computed bv FTan et a,l\ l(200(il) . If the cluster formed in global collapse, 
all th e stars should be to the left of the line. In this case and all the others analyzed by 
iTan et al... (,2006') . the observed age spreads and formation times of clusters are strongly 
incompatible with the idea that star-forming clumps are in a state of global collapse, and 
hence are strongly incompatible with competitive accretion. 

One can also examine the star f ormation rate and compa re to population s of dense gas 
clumps , as i n the r ecent anajysis bv lKrumholz fc Tanll|2006t) . As an example. ICao fc SolomonI 
(l2004fcl) and lWu ~e t al. {20(^ have shown that there is a strong correlation between lumi- 
nosity in the HCN(l-O) line and infrared luminosity that extends from individual cluster- 
forming gas clumps in the Milky Way up to entire ultr aluminous infrared galaxies. Since 
the infrared emission traces the star formation rate (e.g.fkennicutt 1998), and HCN(l-O) 
emission traces the gas mass at densities ti 6 x 10^ H cm" -3 llCao SolomonH2004rfl) . 
a typical density for protocluster gas, this correlation is a direct measure of the star 
formation rate in protocluster gas. lWu et al. (2005), based on observationally-calibrated 
conversions from IR light to star formation rate and HCN(l-O) emission to gas mass 
find that their correlation corresponds to a star formation law « Mhcn/(80 Myr), 
where Mhcn is the mass of gas emitting IICN(l-O). Including some corrections to those 



M. R. Krumholz: Star Formation by Core Collapse 



13 



ge t (10» yr) 



Figure 7. Number of stars out of a sample of 162 versus age in IC348 (histogra m), and dynam- 
ical time in the cluster (dashed vertical line). The histogram is reprinted from lPalla fc Stahlerl 
i2000) with permission from the authors. 



calibrations. iTan et all l)2006(l estimate a star formation law 



(4.5) 



17 Myr 

We ca n compare this observational law to the predictions of theoretical models. FBonnell et all 
l)2003|) simulate a gas clump with a mean density of n = 5.5 x 10* H cm^'^, almost ex- 
actly the same as the density of observed HCN clumps. They find that clusters form in 
a free-fall collapse in which 58% of the gas is incorporated into stars after 2.6 free-fall 
times. At the HCN density of 6 x 10'* H cm~^, the free-fall time is 0.18 Myr, so the star 
formation law one would infer from the simulations is 

inconsistent with the observed relation. In contrast, if one assumes that HCN-emitting 
gas clumps arc virializcd, ctvir ~ 1-5, and turbulent, Mach number M ~ 25, then the 
iKrumholz fc McKea (,2005^) estimate for the star formation rate predicts 

"-WW. 

This is in good agreement with observation, and suggests that HCN-emitting gas can- 
not be in a state of global collapse, either in the Milky Way or in other galaxies. 
IKrumholz & Tan ( 2006) find that repeating this exercise for other populations of star- 
forming objects gives similar results, as shown in Figure |H1 

In addition to the observational data, simulations of star cluster formation including 
feedback are beginning to appear in the litera ture, and these also ca st doubt on the 
idea of global collapse or competitive accretion. iLi fc Nakamural l)2006|) simulate a star- 
forming clump with a code that includes an approximate treatment of the energy and 
momentum injected by outflows from newly formed stars. In their simulations, there are 
no long-lived stagnation points of the type required for competitive accretion. Any such 
stagnation points are rapidly disrupted either by internal feedback from stars within them 
or by external shocking. Energy injected by feedback replenishes energy lost in radiative 
shocks, so the cloud remains virialized and does not collapse. Ins tead, it approaches a 
roughly constant velocity dispersion. In contrast to the results of lBonnell et all l|200,'^ 
without feedback, iLi & Nakamura. (,2006,1 find that only ~ 6% of the mass collapses into 
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Figure 8. Ratio of the free-fall time to the depletion time (defined as the time required to 
convert all gas into stars at the current star formation rate), estimated for gas at a given charac- 
teristic density. The plot shows observationally-determined points {solid squares) f or gas traced 
by CO (l-O) emission (giant molecular clouds), infrared dark clouds taken from the [Simon et al\ 
J2006D catalog, gas traced by H CN(l-O) emiss ion ( Gao & Solomon 20043 IWu et aZ.ll2005D . and 
gas traced by CS(5-4) emission ijShirlev et a^'loOS'). The CS(5-4) point is an upper limit. The 
plot also shows poi nts from simulations (op en di amon ds) bv iKle sscn fc Burkert (2000) and 
IBonnell et al\ |2003) , and the Krumh olz fc McKedl ll2005f) model {hatched region) applied under 

1 - 2 and ID 



the assumption that star formation occurs in clouds wit h virial parameters a.. 
Mach numbers M = 20 ~ 40. The plot is adapted from IKrumholz fc Tanl J200' 



stars per free-fall time, in reasonable agreement with both observational estimates of the 
cluster formation time scale and the star formation rate in dense gas. 

In summary, the question of whether competitive accretion occurs reduces to the ques- 
tion of whether cluster-forming gas clumps are in a state of global collapse, and both 
observational estimates and simulations including feedback seem to rule out this possi- 
bility. The mass in a protostellar core is all the mass that a star will ever have. 



5. Feedback and Accretion 

As shown in the previous sections, massive cores do not fragment strongly, nor do 
the stars that form within them gain significant mass from outside their parent core. 
There remains, however, the problem of getting mass from the core onto a star. This is 
potentially difficult because massive protostars have short Kelvin-Helmholtz times that 
enable them to reach th e main sequence while they are still forming from their parent 
clouds JShu et oi.lll987i) . Once nuclear burning begins, the star will have the immense 
luminosity of a main sequence star of the same mass. Radiation might inhibit accretion 
in two different ways: first, it exerts a direct radiation pressure force on dust grains 
suspended in the incoming gas, and this force could be stronger than gravity. Second, 
ionizing photons can create an HII region around the star, within which the temperature 
is 10"* K and the sound speed is 10 km s^^. The HII region is at very high pressure, 
and may therefore be able to expand and unbind the accreting gas. If either of these 
feedback processes is capable of halting accretion, then the mapping from core mass to 
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star mass will not be simple at the high mass end, and it will not be possible to explain 
the properties of stars in terms of the properties of cores. 

5.1. Radiation Pressure Feedback 

Early spherically symmetric calculations of the effects of radiation pressure on accre- 
tion flows onto massive stars found that ra diation becomes stronger than gravity, and 
hahs accretion, at masses of 20 - 40 Mq ()Kahnlll974 IWolfire fc Cassinellilll987|) for 
typical Galactic metallicities. More recent work has relaxed this limit by considering the 
effects of accretion disks, which allow most of the mass to a ccrete from within the opti- 
cally thick disk, where it is shielded from radiation pressure l|Nakanolll989HNakano et al\ 
Qj)95; Jiiina & Adams 1996). Furthermore, rotational flattening can cause the radiation 
field itself to become asymmetric, as more radiation escapes in the lower density po- 
lar direction than through the higher density equatorial plane. Cylindrically-symmetric 
radiation-hydrodynamic simulations show that this flashlight effect allows accretion to 
slightly more than 40 Mq (|Yorke fc Sonnhalteill200S|) . These simulations also reveal that 
forming a massive star is easier in simulations that include a multifrequency treatment 
of the radiation field than in those that use a gray approximation, probably because the 
coUimation of the radiation field is reduced by the gray approximation. 

More recent 3D simulations show that the flashlight effect persists into 3D, and is ac- 
tually enhanced by a qualitatively new effect. Figure |51 shows time slices from a radiation- 
hydrodynamic simulation of the collaps e of a 100 MfT-,, . 1 pc-rad ius core with an initial 
density profile p oc r~^, performed by iKrumholz et al\ ||2005qI and 2006, m prepara- 
tion). The initial core has a temperature of 40 K and no turbulence, just overall solid 
body rotation at a rate such that the rotational kinetic energy is 2% of the gravitational 
otential energy, a typical rotation rate observed for low mass cores I Goodman et a,I,\ 



1^ 



1993). The simulation shows that, at stellar masses ^17 Mq, stellar radiation is too 
weak to reverse inflow, and the gas falls onto an accretion disk and accretes in a steady, 
cylindrically-symmetric flow. At about 17 Mq, stellar radiation begins to reverse inflow 
and drive radiation bubbles above and below the accretion disk. The interiors of these 
bubbles are optically thin and very low density, while the walls reach densities ^ 10^° 
H cm~'^. Gas reaching the bubbles ceases to move radially inward. Instead, it travels 
along the bubble wall until it falls onto the accretion disk. At that point, it is shield from 
stellar radiation and is able to accrete onto the star. As shown in Figure [Till the bubbles 
also collimate radiation, increasing the flux in the polar direction and decreasing it in 
the equatorial plane. At the time shown in Figure [TUk . the radiation flux just outside the 
radiation bubble in the polar direction exceeds that in the equatorial plane at the same 
radius by more than an order of magnitude. 

The bubbles expand asymmetrically with respect to both the polar axis and the equa- 
torial plane due to an instability. The direction in which a bubble has expanded farthest 
and cleared out a low optical depth region represents the path of least resistance for 
radiation leaving the star to escape to infinity. As a result, radiation is coUimated in 
whichever direction the bubble has expanded farthest. In turn, this causes the radiation 
pressure force to be largest in that direction, which amplifies the rate of bubble expan- 
sion. Figure^ shows this asymmetric expansion phase. As the simulation progresses the 
instability becomes more violent and the bubbles start to collapse, as shown in Figure 
The collapse is due to radiation Rayleigh- Taylor instability: the gas is a heavy fiuid, 
and the radiation is a light fluid trying to hold it up, an unstable configuration. Once 
the bubbles collapse, as shown in Figure Et, inf ailing gas is deflected by radiation onto 
the dense remnants of the bubble walls, which are optically very thick and therefore 
self-shielding against radiation. The gas travels along the walls onto the disk and then 
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Figure 9. Slices in the XZ plane showing t he density (grays c ale) an d velocity (arrows) in a 
simulation of the collapse of a 100 Mq core by lKrumholz et aJU (12005 of) . The panels show times 
of (a) 1.5 X lO** yr, (b), 1.65 x IC* yr, and (c) 2.0 x 10* yr after the start of the simulation. 
The corresponding stellar masses at those times are 21.3 Mq, 22.4 Mq, and 25.7 Mq. Note that 
these masses are considerably larger than those of the turbulent runs discussed in § |3] at the 
same times primarily due to the steeper initial density profile assumed in these simulations. 
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Figure 10. Slices in the XZ plane showing the density (grayscale) and radiation flux multiplied 
by 47rr^ (arrows) in the simulation shown in Figure|5| For clarity, the flux vectors in the optically 
thin bubble interior in panel (a) have been omitted. The times shown in panels (a) and (b) are 
the same as those shown in panels (b) and (c) of Figure]^ 

accretes. In this configuration the remnant walls collimate the accretion, as shown in 
Figure [TUb . The radiation flux behind the walls is very small, so gas accreting from that 
direction feels very little radiation pressure force. 

Simulations to date have reached « 34 Mq in this manner, with 5—10 Mq more in 
the accretion disk and no sign of a r eversal of i nfall. Note that these calculations use 
gray radiative transfer, for which case lYorke fc S onnhaltcr (2002) found an upper limit 
of about 20 Mq for the stellar mass. This suggests that multifrequency 3D calculations 
would show even easier accretion. 

iKru mholz et al. (2005c) point out a second mechanism for overcoming radiation pres- 
sure feedback. Massive protostars have outflows that appear to be siniply sc aled-up ver- 
sions of the outflows from lower mass stars (see iBeuther fc Shet>herdll2005l for a recent 
review), or that are somewhat wider in angle for the most massive, O type stars. These 
outflows can have a very strong effect on the radiation field. Outflows are launched from 
the inner accretion disk around the star, where the gas is so hot that any dust grains 
present within it sublime. As a result, outflows contain almost no dust when they are 
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launched, and simple calculations show that due to the velocities of ^ 500 km s^^ with 
which they flow, there is insufficient time for large dust grains to reform before the out- 
flowing gas reaches distances ^ 0.1 pc from the source star. This makes the outflow 
cavities optically thin, so they can coUi mate radiation and carr y it out of the optically 
thick accreting envelope very effectively. iKrumholz et al\ ()2005cf) use Monte Carlo radia- 
tive transfer calculations to explore this effect, and show that it can reduce the radiation 
pressure force on infalling gas by an order of magnitude. This can in turn shift the ra- 
diation pressure force from a regime where it is stronger than the gravitational force to 
one where it is weaker than gravity, allowing accretion to continue where it might have 
been halted without the protostellar outflow. 

It is unclear how this effect will interact with radiation pressure bubbles, since the 
simulations in which the bubbles appear do not yet include a model for protostellar 
outflows. Nonetheless, the overall conclusion one may draw from the two effects is that 
radiation pressure cannot easily halt accretion. There are at least two mechanisms which 
are capable of rendering the radiation field aspherical, allowing accretion over the range 
of solid angle where the radiation pressure force is reduced. Gravity is very effective at 
channeling gas into those solid angles, and then onto the star. 



5.2. Ionization 

The second form of radiative feedback that accretion onto a massive star must over- 
come is ionization, firs t recognized as a problem by iLarson fc Starrfieldl l)l97l|) and 
lYorke fc Kru'^ ^97^. If the ionizing photons from a massive protostar are able to 
reach gas at distance from the star where the escape speed is <^ 10 km s~^ then the ion- 
ized gas will be unbound and may escape to infinity. However, IWalmslevI l|l995l) shows 
that sufficiently rapid accretion can avoid this problem by keeping the ionized region 
quenched against the stellar surface. For spherically symmetric infall, if the accretion 
rate is larger than 



where S is the ionizing luminosity of the star in photons and a(2) is the recombination 
coefficient to excited levels of hydrogen, then all ionizing photons will be absorbed near 
the stellar surface. The escape velocity there is much larger than 10 km s~^, so this 
ionized gas will be confined by stellar gravity and accretion may continue unimpeded. 

Observational estimates of HII region lifetimes indeed appear to require confinement of 
this sort. Ultracompact HII regions, roughly those <J 0.1 pc in radius, have a dynamical 
expansion time 0.1 pc/10 km s^^ ~ 10"* yr, but the number of o bserved ultracompact 
HII regions implies tha t their lifetimes must be closer to 10^ yr ijWood fc Churchwelll 
llflSatlKurtz fi^ oi.lll994l) . Moreover, some observed syst ems show sign s of both accretion 
and t he presence of an ultracompact HII region fe.g. iHofner eta!\ K999: So Uins et all 
l2005tl. providing dire ct evidence that accretion can continue past HII regior i forrnation. 

iHoare et al\ |2005|) gives a recent review of ultracompact HII regions, and lKetol l)2006() 
attempts to construct a unified theoretical framework within whi ch one can understan d 
the relationship between HII regions and accretion (though see lBeuther et a/.l 12005 oj) . 
In this model, the behavior of an HII region is controlled by the relative sizes of the 
ionization radius within which the ionization and recombination rates balance and all 
ionizing photons are absorbed, the gravitational radius where a spherical accretion flow 
accelerates from subsonic to supersonic relative to the sound speed in ionized gas, and 
the disk radius, at which angular momentum in accreting gas causes infall to become 
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aspherical and a disk to form. This model can roughly explain observed HII region 
morphologies. One of its consequences is that if the ionization radius is smaller than 
the gravitational radius, the HII region will remain confined and accretion will continue 
unimpeded. For the accretion rates of lO"** — 10~^ Mq yr~^ one expects given the 
conditions in massive cores, massive protostars will remain in this regime until they have 
exhausted essentially all of their available gas. Thus, ionization feedback does not present 
a significant barrier to accretion. 



6. Conclusions and Future Directions 

6.1. Massive Star Formation from Massive Cores 

The results summarized in the previous sections combine to give a unified scenario of 

massive star formation: massive stars form by the gravitational collapse of massive cores. 
We observe these cores, from a few tenths to a few hundred Mq in mass, in the dense 
molecular clumps where star clusters form. Their mass distribution matches the stellar 
IMF, and their spatial distribution is roughly consistent with the spatial distribution of 
massive stars in clusters. When massive cores collapse, they do not fragment strongly 
because radiation feedback from the high accretion rates they produce warms their inner 
regions, raising the .leans mass and suppressing fragmentation. As a rcsiilt, all of their 
mass falls onto one a few stars. This occurs even though the cores are turbulent and 
contain complex density and velocity structures. 

The stars that form in massive cores do not gain significant mass from outside their 
parent cores, because gas that is not gravitationally bound to them at birth is too tur- 
bulent for significant accretion to occur. The mostly likely agent for keeping the gas 
turbulent is feedback from protostellar outflows, and both observational estimates and 
simulations suggest that this mechanism works. Regardless of the agent, though, obser- 
vations of the age spread of young clusters and the star formation rate in pre-cluster gas 
unambiguously require that star cluster form from gas that is not in free-fall collapse. 

The gas falling onto a massive protostar at the center of a massive core must accrete in 
the face of significant radiation feedback in the form of radiation pressure on dust grains 
suspended in the inflowing gas, and ionizing radiation that generates high pressure ion- 
ized gas. However, neither mechanism substantially impedes accretion. Radiation escapes 
through low optical depth channels, which either form spontaneously due to radiation- 
hydrodynamic instabilities or are created by protostellar outflows. Gas accretes through 
high optical depth channels onto an accretion disk, where it is shielded from protostellar 
radiation. At the same time, the accretion flow is capable of absorbing all the ionizing 
stellar photons near the star, where the escape speed is larger than the ionized gas sound 
speed. As a result, formation of an HII region is either suppressed entirely, or the ionized 
region is kept close to the star where gravity confines it and prevents it from expanding. 
Accretion continues through the ionized region. 

Because gas from a massive core accretes onto one or two stars, with no significant 
impedance from feedback, and those stars do not gain a significant amount of mass 
from outside their parent cores, there is a direct mapping from the properties of cores 
to the properties of stars. Consequently the mass and spatial distributions of young 
clusters are direct imprints of the mass and spatial distributions of cores in their gas 
phase progenitors, possibly with some additional redistribution of mass due to dynamical 
mass segregation over the several crossing times that it takes for a cluster to assemble. 
Observables in the gas phase connect directly to observables in the stellar phase. 
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6.2. Problems for the Future 

While this scenario explains a great deal of what we observe about massive stars, it also 
omits a number of physical effects that a more complete picture should include, and it 
needs more rigorous testing against observation. In this section, I discuss the prospects 
for improvement in these areas in the next few years. 

6.2.1. Magnetic Fields 

None of the models of massive star formation proposed to date have included the ef- 
fects of magnetic fields on the dynamics of the collapsing gas. Partly this is due to a 
lack of observational information about magnetic fields in massive star-forming regions. 
ICrutchen (,.2QQ5) summarizes the state of observations, and concludes that magnetic fields 
are marginally dynamically significant in massive star-forming cores. However, this de- 
termination is extremely uncertain due to the difficulty of interpreting observational 
indicators of magnetic field strength. Part of this uncertainty comes from geometry. The 
critical magnetic flux required to hold up a core against gravity depends on the core 
shape and mass distribution, and differing assumptions about the shape can produce 
qualitatively different conclusions about how th e magnetic field str ength compares to 
that required for it to be dynamically significant ijBourke et aZjEoOlj) . 

Beyond the geometric effect, measurements of field strengths themselves are quite un- 
certain. The two most commonly-used techniques for determining magnetic field strengths 
in molecular gas are Zeeman splitting of OH or CN lines, and the Chandrasekhar- Fermi 
method in which one uses the dispersion of polarization vectors seen in dust emission to 
estimate how the ratio of magnetic to kinetic energy density. Both of these techniques 
suffer from major systematic errors. Zeeman splitting averages the magnetic field along 
the line of sight, potentially washing out the signal, and is also only sensitive to regions 
where the observed species are present. Since m any species freeze out onto dust grains in 
the densest parts of cores (|Tafalla et allhOO'^ . Zeeman measurements may reveal more 
about magnetic field strengths in the diffuse gas around cores than in cores themselves. 
For Chandrasekhar- Fermi measurements, one also faces an uncertainty about where along 
the line of sight the observed polarized light is being emitted, and whether one is really 
measuring field in the core or in its outskirts. This uncertainty is made even worse by 
the fact that the magnetic field one deduces for a given dispersion of polarization vectors 
depends on the gas density, so uncertainty in the gas properties in the emitting region 
translates directly into uncertainty about the magnetic field strength. 

The level of uncertainty is driven home by the result that, in at least some regions where 
both have been used, the inferred ma gnetic field strengths differ by a factor of ~ 100 (e.g. 
ICrutchedllQQfltlLai et qiJl2nOlL 120021) . While in principle this is possible because Zeeman 
splitting measures the field along the line of sight and Chandrasekhar-Fermi measures 
the field in the plane of the sky, it seems implausible that the magnetic field in multiple 
sources could be so perfectly aligned perpendicular to the line of sight. A more likely 
explanation is that there are large systematic errors in one or both techniques, or that 
the two techniques are simply measuring the field in different regions. 

In the absence of more conclusive observational data, it is hard to know whether it is 
important to include magnetic fields in massive star formation models or not. In prin- 
ciple one could make models both with and without dynamically significant magnetic 
fields to explore their effects, but at present magnetohydrodynamics codes using either 
smoothed particle hydrodynamics or adaptive mesh refinement, one of which is required 
to achieve the requisite dynarnic range to simulate star f ormation, are still in their in- 
fancy fe.g. IPrice fc Monag-haiJEoO^ ICrockett et aJlhOO^ . There are no codes past the 
experimental stage that include both MHD and radiative transfer, which the discussion 
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above demonstrates is absolutely crucial. The question of the influence of magnetic fields 
will therefore have to wait for progress in both code development and observations. 

6.2.2. Improved Simulation Physics 

In addition to magnetic fields, there are numerous other physical effects that the cur- 
rent generation of simulations ignore. It is chastening to realize that almost every new 
piece of physics that has been added to simulations of massive star formation has revealed 
a qualitatively new and unexpected phenomenon: adding radiative transfer to simulations 
of turbulent cores shows that fragmentation is much weaker than one might suppose from 
purely hydrodynamic calculations (§01; moving from gray to multifrequency radiative 
transfer in two dimensional simulations doubles the maximum stellar mass formed due 



three dimensions with gray ra diative transfer r eveals the formation of unstable, asymmet- 
ric radiation bubbles (,Krumholz et "aIl l2Q05Ql and § El; includ ing protostellar outflows 
in mo dels greatly reduces th e efficacy of radiation feedback on the scale of individual 
cores ijKrumholz et a/.ll2005cj) . and changes the dynamics of star cluster formation from 
rapid collapse and competitive accretion to quasi-equilibriu m behavior where tur bulent 
decay is balanced by injection of kinetic energy by outflows l)Li fc Nakamurall2006|) . This 
pattern suggests that there may be new effects yet to be discovered as simulation physics 
improves. An important corollary of this is that we have very likely reached the limit of 
what pure gravity plus hydrodynamics simulations with no feedback can teach us. Major 
work in the future will not come from simply running more and bigger simulations of the 
same type, but from adding new physics to the problem. 

The next logical steps in modeling massive star formation include extending multi- 
frequency radiative transfer to three dimensions, using radiation hydrodynamic codes 
to simulate the formation of star clusters instead of just individual stars, and includ- 
ing protostellar outflows in radiation-hydrodynamic simulations of massive cores. These 
are all relatively straightforward extensions of existing techniques. The primary limit to 
using them is computer power, since gray radiation-hydrodynamic simulations in three 
dimensions take months on available supercomputers. However, Moore's Law will help 
in this regard, and these simulations should become feasible in the next few years. More 
difficult are the problems for which not only computer power but also our current sim- 
ulation methods are inadequate, and for which we will have to develop new techniques. 
One problem in this category is putting more realistic dust physics into three dimen- 
sional models, including coagulation, shattering , and differential motion of grains and 
gas, which may substant ially change opacities l|Preibisch et all llQQSt ISonnhalter et all 
Suttner et al Another is radiative transfer that goes beyond flux-limited 

diffusion (e.s iHaves & Norman 2003). Improvements in these areas are likely to take 
somewhat longer, perhaps becoming reasonable on the time scale of five or more years. 

6.2.3. Connecting Models to Observations 

A third area in which substantial progress must be made is connecting simulation and 
theoretical models to observations. There has been a fair amount of work on relating 
simulation results to statistical diagnostics of young clusters such as the IMF and binary 
and brown dwarf properties. However, these indicators are quite difficult to use as a 
means of distinguishing theories because they focus on the outcome only after the star 
cluster is fully formed, and there appear to be mu ltipl e mechanisms capab le of producing 
the same ou tcome fe.g.lPadoan fc Nordlundl2002l andlBonnell et 'aZll2003l on the origin of 
the IMF, or iRonnell R a,tfl l2nn,'j a,nd iKnimholz et oiJl2006r;l on the binarity properties 
of massive stars). In addition, observations of the results of star formation are difficult to 
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connect to models because we do not know the properties of the progenitor gas system 
for a given steUar system. It is often possible to reproduce a given stellar system simply 
by tuning the properties of the assumed progenitor. 

Better diagnostics are likely to come from comparing the results of simulations and 
theoretical models to observations of systems that are still actively forming stars. Masers 



provide one potential connection, sin ce they make it possib le to trace gas properties and 
kinematics on very s mall scales (e.g . 'Torre lles et a^JlioOll who find expanding bubbles 
traced by masers, and lGreenhill et al. 2004. who detect a disk and an outflow cavity traced 
by masers). These may be particularly useful for tracing radiation-driven bubbles, since 
the density of 10^ — 10^'^ cm~^ and temperature of a few hundred K of bubble walls is fa- 
vorable for maser emission. As discussed elsewhere in this volume, observational searches 
for accretion disks around massive stars, which are expected for massive star formation 
by core collapse but probably not for coUisional formation, are another exampl e. Mor- 
pholo gies of warm gas, as seen in high resolution observations of "hot cores" fe.g. lGaravl 
provide another potential point of comparison to simulations. One can carry such 
comparisons even further by post-processing the simulation results with radiative trans- 
fer codes to produce simulated line profiles and intensity maps. On e can also look for 
tellta le signs of star formation by collisions, such as infrared fiares ijBallv Sz Zinnecke^ 
l2005l) or embedded clusters in which the stellar density reaches the ~ lO"^ pc ^ required 
for efficient collisions - such high densities should produce d istinctive spectral energy dis- 
tributions that peak at wavelengths approaching 100 /zm l|Chakrabarti &: McKeell2005L 
and Chakrabarti & McKee, 2006, in preparation). Yet another technique is to use extra- 
galactic observations, in which one can o bserve the entire galacti c population of a given 
sort of object (e.g. HCN-emitting clouds iG^^^^^IsfflS^Jl^^^L ^^^"^ then make argu- 
ments based on population statistics (e.g. iKrumholz fc Tanll2006l) . Any or all off these 
are likely to be far more definitive than comparisons to stellar systems. 

Prospects for improvement in this area of massive star formation theory are quite 
good. Current and planned millimeter instruments such as PdBI, the SMA, CARMA, 
and ALMA will allow higher resolution examination of more distant, embedded star- 
forming systems than has heretofore been possible, and surveys or both Galactic and 
extra-Galactic star-forming regions with Spitzer are providing large databases that we 
can use to test models on the level of populations. These sources of data have already 
begun to let us distinguish between models of massive star formation, and will provide 
even more powerful diagnostics in the next few years. 
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